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Restrictions on the injection energy of positrons 
annihilating near the Galactic center 
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ABSTRACT 

The origin and properties of the source of positrons annihilating in the Galactic Center 
is still a mystery. One of the criterion, which may discriminate between different 
mechanisms of positron production there, is the positron energy injection. Beacom 
and Yiiksel (2006) suggested a method to estimate this energy from the ratio of the 
511 keV Une to the MeV in-flight annihilation fluxes. From the COMPTEL data they 
derived that the maximum injection energy of positron should be about several MeV 
that cut down significantly a class of models of positron origin in the GC assuming 
that positrons lose their energy by Coulomb collisions only. However, observations 
show that the strength of magnetic field in the GC is much higher than in other 
parts of the Galaxy, and it may range there from 100 jiG to several mG. In these 
conditions, synchrotron losses of positrons are significant that extends the range of 
acceptable values of positron injection energy. We show that if positrons injection in 
the GC is non-stationary and magnetic field is higher than 0.4 mC both radio and 
gamma-ray restrictions permit their energy to be higher than several GeV. 
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1 INTRODUCTION 

One of the interesting and still unsolved problems is the 
origin of 511 keV annihilation emission from the Galac- 
tic Bulge. It is observed as an extended diffuse emission 
from 5° — 8° radius region with the flux ~ 8 x 10~* ph 
cm~^s~^ tha t requires the rate of po s itron production there 
~ 10^^ s~^ (iKnodlseder et al.l l2005l : IChurazov etHI l2005l : 
I Jean et "all 120061 ) . These observations showed that the en- 
ergy of annihilating positrons was about 1 eV. On the other 
hand, all potential sources of po sitrons in the Galaxy like 
SN stars (|Knodlseder et al.ll2005j V massive stars generating 
the radioactive '^^Al ( Prantzoi |2006|) . secondary positrons 



High energy positrons annihilate "in-flight", thus pro- 
ducing continuum emission in the range E > 511 keV. 
A prominent 511 keV line emission is generated by these 
positrons when their energy is decreased to the thermal one 
due to energy losses. Therefore, one can expect that the 
continuum and line emission are proportional to each other. 

For the lifetime of in-flight annihilation, Tif, and the av- 
erage cooling time Tci of high energy positrons, one can esti- 



mate the expected flux of in-flight annihilation, 
from the observed 511 keV line emission, lune, as 



from p-p collision s ( Cheng et al.l 20071 . le pton jets of 



-'^i-so MeV 
where 



if Tif » Tel 



1-30 MeV 



(1) 



AGNs (lTotanill2006l'). dark matter annihilation (|Boehm et all 
l2004l: [S izun et al. I l2006l) mic r oquasars and X-ray b i naries 
( Weidenspointner et all l2008l : iBandvopadhvav et all bOOTl ) 



etc. generate positrons with energies > 1 MeV. This means 
that positrons should effectively lose their energy before an- 
nihilation and thus generate emission in other than 511 keV 
energy ranges. Therefore, the injection energy of positrons is 
an essential parameter for modeling annihilation processes, 
and it can be in principle discriminated from observations. 



nt = [na,f{E)v{E)]-^ 
and 



Tel = 



dE 



(dE/dt)^ 



(2) 



(3) 



Here n is the plasma density, and Ei^j, Eth and v{E) are the 
injection energy of positrons, the energy of thermal plasma. 
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and the positron velocity, respectively, a"i/ is the cross- 
section of in-flight annihilation. The function (dE/dt)ci is 
the rate of energy losses defined as sum of Coulomb, syn- 
chrotron, inverse Compton, bremsstrahlung etc. losses: 



{dE/dt)^, = (dE/dt)^^^^ + (dE/dt)^^^ 
+ {dE/dt),c + {dE/dt)^,. + ... 



(4) 



If cooling of the positrons is only due to the Coulomb 
losses, then r^/ and Td are proportional to n~^, and the rela- 
tion ([T]) is independent of the medium density. So, this ratio 
of the continuum in-flight and the annihilation emission is 
universal and ca n be applied even to a me dium with an 
unknown density. iBeacom and Yiiksell (|2006l ) assumed that 
positrons in the Galactic center (GC) loose their energy by 
Coulomb interactions only, and they suggested to use this 
ratio for the analysis of the annihilating positron origin in 
the GC. In the above-mentioned models the injection en- 
ergy of positrons is expected in the range from several to 
hundreds MeV. Therefore, the in-flight gamma-ray emission 
is also expected in this energy range. 

The MeV flux from the cen tral part of the G alaxy 
was observed by COMPTEL fsee. lStrong et al.lll998l ). The 
origin of this emission is still unclear since the known 
processes of gamma-ray production (like inverse Compton, 
brem sstrahlung etc.) are unable to g enerate the observed 
flux (IStrong et all l2005l : iPorter et al. 2008 ). ICheng et all 
l|2007^ assumed that this excess in the GC direction might be 
due to the in-flight annihilation of fast positrons. However, 
it is observed not only in the direction of the GC. The ex- 
cess i s almost constant along the Galactic disk (jStrong et al.l 
1 19981 ) where the intensity of annihilation emission is lower 
than in the Galactic centre. This makes problematic the in- 
flight interpretation of this excess in the disk since the ratio 
511 keV flux/in-flight continuum is constant. If the in-flight 
flux is responsible for the MeV excess in a relatively narrow 
central region (< 5°), absolutely the sa me excess in other 
parts of the disk remain s unexplained (see lSizu"ii et al.ll2006l : 
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Therefore, in 



20081) 



Beacom and Yiiksell (|2006D and latter in 



ISizun et all l|2006h a more firm constraint on the in-fiight 
gamma-ray flux from the Galactic center was suggested. Ac- 
cording to their criterion the in-flight flux should not exceed 
several statistical errors of the COMPTEL measurements. 
That gives an upper limit for the injection energy about 
several MeV. Then models assuming higher injection energy 
should undoubtedly be rejected. 

Below we show that under some conditions the injection 
energy may be higher than 10 MeV in contrast to conclusions 
made in papers mentioned above, and, thus, there is a room 
for models assuming inj e ction o f high energy positrons. 

Thus, ICheng et all l|2006l . Ho^ assumed that these 
positrons are secondary and generated by collisions of rela- 
tivistic pr otons injected from bl ack hole jets. The theoretical 
analysis of llstomin fc Soli (|2009D confirmed the hadronic ori- 
gin of jets and showed that protons were accelerated there by 
the stochastic and the centrifugal acceleration up to energies 
E„ ~ 10^" eV that might offer an explanati on to the recent 
resul ts of the Pierre Auger collaboration (jAbraham et al.l 
l2007l ). If such or similar mechanism produces indeed enough 
relativistic protons with Lorentz factor 7 > 2 in the vicinity 
of the central black hole then we do expect there an effec- 



tive production of secondar y positrons w i th energies above 
30 MeV, just as assumed in lCheng etahl l|2006l . |2007| ). 

Processes of p — p collisions produce also a flux of 
gamma-rays in the range above 100 MeV by decay of 
7r°-meson, and below this energy by, so-called, internal 
bre msstrahlung radi ation of secondary electrons (see for de- 
tail lHavakawalTl964l ). A flux of gamma-rays in the 1 to 30 
MeV range from internal bremsstrahlung may be higher 
than the mentioned in-flight flux. Thus, iBeacom et al.l 
(|2005l ) showed that the internal bremsstrahlung flux is very 
significant, if positrons in the GC are generated by dark mat- 
ter annihilation. However, in the dark matter model positron 
production in the GC is stationary. On the other hand, from 
the restrictions derived from EGRET data it follows that 
the positron production in the GC should be strongly non- 
stationary, if these positrons are generated hy p — p collisions 
( Cheng et al. 2006, 2007). The fiux of gamma-rays from p—p 
collisions is significant during a very short period after a star 
accretion onto the black hole. At present this fiux has de- 
creased in several orders of magnitude from its initial value 
and, therefore, is unseen. 

Below in section [3] we shall show that the condition of 
non-stationarity is also required to fit radio observations. 



2 MEDIUM PROPERTIES IN THE VICINITY 
OF THE GALACTIC CENTER 

As follows from observations, the central 200 pc region of 
the Galaxy is strongly nonuniform. The inner bulge (200- 
300 pc) contains (7 — 9) x 10^ Mq of hydrogen gas. In spite 
of relatively small radius this region contains about 10% 
of the Galaxy's molecular mass. Most of the molecular gas 
is contained in very compact clouds of mass 10* - 10*^ Mq, 
average densities ^ lO^cm""^. 

However, this molecular gas occupies a rather small part 
of the cent ral region, most of wh ich is filled with a very hot 
gas. ASCA lKovama et al.l l|l996f ) measured the X-ray spec- 
trum in the inner 150 pc region which exhibited a number of 
emission lines from highly ionized elements which are char- 
acteristics for a 8— 10 keV plasma with the density .4 cm"''. 
Later on Chandra observations o f lMuno et al] l|2004 ) showed 
an intensive X-ray emission at the energy Ex ~ 8 keV from 
the inner 20 pc of the Galaxy. The plasma density was es- 
timated in limits 0.1 — 0.2 cm~^. Recent SUZAKU mea- 



surements of the 6.9/6.7 keV iron line ratio ama et al.l 
I2OO7I ) was naturally explained by a thermal emission of 6.5 
keV-temperature plasma. 

One should note that there is no consensus on the mag- 
netic field strength in the GC. Estima tions ranges from 
about or smaller than hundred mG (see ISpergel and Blit j 
1 19921 : iLaRosa et ai] 12005 ^: 'Higdon et al '200^) , up to sev- 
eral m C ( Plante, Lo and Crutchor 1995; Yuscf-Zade h et al.l 
1 19991 . see in this respect the review of lFerrierell2009l ). Radio 
observations of the central regions show that the structure of 
strong magnetic fields is nonuniform and it concentrates in 
filaments which extends up to 200 pc from the GC. The 
region containing mG magnetic field is estimated by the 
angular size 1.5° x 0.5° (see e.g. Morris and Serabynl 19961 : 



Lbyr 



lYusef-Zadeh. Hewitt and Cottonll2004l : lMorrij|2006 

In these magnetic fields synchrotron losses are essential 
even for positrons with energies 3-30 MeV. One can see from 
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Figure 1. Characteristic cooling times of positrons for different 
values of magnetic field strength. 



Fig.[l]that if the magnetic field strength is as high as 3 mG, 
the cooling time of 1 GeV positrons is the same as that 
of positrons with energies about 1 MeV experiencing only 
Coulomb losses. 

Structure of the magnetic field and the way electrons 
and positrons interact with it is also unknown. In our model 
we assumed that positrons interact with magnetic field 
structures violently thus allowing us to average the magnetic 
field over central 100 pc radius. However if the interaction is 
we ak and can be neglec t ed th e original constraints deduced 
bv lBeacom and Yiiksell (|2006! ) prevail. 

We shall show that the injection energy of positrons can 
be much higher than 1 MeV because of synchrotron losses 
in the Galactic center that extends significantly the class of 
models explaining the origin of annihilation emission from 
the bulge. 

If it is not specified we accepted below that the cen- 
tral 100 pc radius region is filled with strong magnetic field 
and the gas density in the central region (r < 500 pc) is 
> 0.2 cm~^. The rate of ionization losses depends on the 
medium ionization degree. Therefore, we consider two cases 
of neutral and fuUy ionized medium. 



3 EFFECTS OF STRONG MAGNETIC FIELDS 
IN THE GC 

If positrons propagate through the region of strong magnetic 
fields in the GC, they effectively lose their energy and emit 
radiation in the radio range of frequencies. Thus, for H = 3 
mG positrons with energies 10-100 MeV generate emission 
in the range 



1^0 



: 0.29 



3eH 



E 



1 - 100 MHz 



(5) 



4nmc \mc^ 

For the magnetic field of 3 mG the rate of synchrotron 
losses of relativistic electrons and positrons is higher than 
for any other process of losses. Therefore, an intensive syn- 
chrotron emission is expected from this region. 

Observations of the GC gave t he following values 
of radio fiux from there (see, e.g., iMezger et al.l 1 19961 : 
iLaRosa etal] 120051 ). This flux is about 10 kjy in the fre- 
quency range 300 - 700 MHz for the central region 3° x 2°. 



For the region 1.5° x 0.5° this flux is almost one order of the 
magnitude smaller at the frequency 330 MHz, i.e. ~ 1 kJy. 

On the other hand, simple estimations of synchrotron 
emission shows that if electrons and positrons are injected 
with energies about 100 MeV, then the radio flux from the 
GC is much higher than 1 kJy. 

Below we demonstrate this point from very simple esti- 
mates of the radio emission from the GC. In the stationary 
state the total spectrum of positron in the GC injected at 
the energy Eo with the rate Q is 

Q 



N{E) 



-e{Eo - E) 



(6) 



\dE/dt\ 

where the rate of synchrotron losses is 

dE_ _ 2 e^H^ ( E 
dt 3 c^m'^ 

Here and below N{E) is the density of electrons and 
positrons with energy E. The radio flux from the GC at 
Earth can be calculated from the integral 



\mc-' J 



(7) 



N{E)P{u,E)dE 



(8) 



where E) is the emissivity of a single positron with the 
energy E, and d = 8 kpc is the distance between Earth and 
the GC. To get an e stimation we use the app roximation for 
the function P from iBerezinskii et all (|l990h 



P{u,E)c^^S{u^uo{E)) 

Then we have the the radio flux at the frequency v 

Q mc^ I — — 
8^ — ^"/^ 



(9) 



(10) 



where v — -^0.29 • SeH/Ainnc. For the production rate of 
secondary electrons Q — 2 ■ 10'*'^s~^, this equation gives a 
flux at Earth about 10^ Jy at the frequency u ~ 100 MHz 
that is several order of magnitude higher than observed. This 
means that either the magnetic field strength and the injec- 
tion energy of positrons should be smaller than those used 
in this estimation or the situation is strongly non-stationary. 

For estimations of the non-stationary case we use the 
source function in the form 



Q{E,t) = NoS{t)S{E - Eo) 



(11) 



where A'o is the total number of injected positrons. Then 
the distribution function of positrons is (see section [5]) 



N{E,t) = NoS[E 



Eo 



PEot + 1 



(12) 



where /3 is the coefficient in Eq. (O, and the radio flux has 
the form 

iVo dE 



v{E ■ 



t) 



dt 



(13) 



As one can see that the current energy of electrons is inde- 
pendent of their injection energy Eq and equals E ~ l//3t 
for long enough t. 

From Eq. (|13p it follows that the radio flux just after 
the injection of positrons equals ~ 3 • 10^ Jy for the A^o — 
2 • 10^^ in each capture events that is necessary to produce 
the observed annihilation flux from the GC. However, for the 
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time t ~ 10^ years (the average period between str captures) 
the peak of intensity shifts from the frequency 100 MHz 
to several MHz where this radio flux cannot be observed 
because of absorption in the interstellar gas, and just this 
effect is a key point of our analysis presented below. 

As an example we show in Fig. [2] the spectrum radio 
emission from the GC region produced by secondary elec- 
trons and positrons at the time 10^ yr after the capture. The 
magnetic field strength in the GC is 2 and 3 mG. As one can 
see for long time after the capture the peak of emission shifts 
from the frequency of hundreds MHz to several MHz, and 
the intensity of radio emission is negligible at 330 MHz. For 
calculatio ns of radio emission we used the accurate equation 
(from e.g. iBerezinskii et al.lll99(il ). 

oo oo 
v/v^ 

where N{E) is the total number of electrons and positrons 
with the energy E, d — S kpc is the distance to the GC, and 
Ka(x) is the McDonald function. 

In Fig. Owe presented the limitations of positron injec- 
tion energy Eq derived from radio data for different values 
of the magnetic field strength and different times T passed 
from the last star capture. We considered two possible cases: 
the injection in the form of delta- function, Eq. (Illf) . shown 
in the figure by solid lines and power-law injection spectrum 
(dashed line): 

No{E) = No5{t)E-^e{E - Eo) . (15) 

One can see that if the injection is non-stationary, then the 
injection energy of positrons can be extremely high in case of 
strong magnetic fields in the GC. Indeed for almost station- 
ary situation when characteristic period between injections 
is small T ^ 1 x 10** yr the maximum allowed energy in case 
of J/ ^ 1 mG cannot exceed 30 MeV due to radio limita- 
tions. In case of power-law spectrum situation is even worse. 
However if period is long enough and magnetic field is strong 
the situation changes: for T = 1 x 10^ yr and H — 2 mG the 
maximum energy is about 1 GeV and it can be even higher 
for longer periods and higher values of the magnetic field. 




Figure 2. Radio flux at Earth produced by positrons from a 
single capture event at the time 10^ yr after the capture. The 
strength of magnetic field H = 3 mG (solid line) and H = 2 mG 
(dashed line) 




10* 10"^ 
H(G) 

Figure 3. Limitations from radio observations for different val- 
ues of time T passed since capture. In ascending order T = 1 X 10'*, 
5 X 10*, 10®, and 2 X 10® years. Solid lines correspond to injection 
in form of delta-function, dashed lines correspond to power-law 
injection. 



4 CONSTRAINTS OF THE MODEL OF 
ANNIHILATION EMISSION 

There are restrictions for our model which follow from radio 
and gamma-observations: 

(i) The fiux of radio emission from the GC at present 
should not exceed t he value of 1 kjy in the range 300 - 
700 MHz (see, e.g., iMezger et al] llOQel ) that restricts the 
number of high energy positrons in the GC. The radius of 
high strength magnetic field sphere was taken to be rn ~ 
100 pc; 

(ii) The flux of annihilation emission observed at Earth 
from this region is 8 • 1 0""^ photons cm~'^s~ ^ for the 
FWHM=6° central region (jChurazov et al.ll2005l ): 

(iii) The fiux in the MeV energy range was measured by 
COM PTEL. As in iBeacom and Yuksell (|2006l '): ISizun et all 
(|2006l ) we accept that the flux of in-fiight annihilation emis- 



sion should not exceed the COMPTEL 2cr level of in the 
energy range 1-30 MeV, if it is due to in-flight annihilation. 

These restrictions are presented in Table [T] 



Table 1. Constraints of the model 



<I.^(i/ = 330 MHz) 


-^^1-30 MeV 


^511 keV 


1.5° X 0.5° region 




FWHM =6° 


1 


< 2cr of COMPTEL 


8 ■ 10-4 


kJy 




ph cm^^s^^ 
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5 EVOLUTION OF THE SPECTRUM OF 
RELATIVISTIC POSITRONS 

We start from the spatially uniform model, which estimates 
the total flux of gamma-rays and radio emission (i.e. in- 
tegrated over the volume of emission). The evolution of 
positron spectrum can be described by the equation 

dN d /dE. 
'dt'^dE [~dt^ 

where dE/dt is the rate of synchrotron losses (see Eq. ([Til). In 
case of instantaneous injection Q{E,t) = No{E)5{t) where 
No{E) is the injection spectrum of positron, the solution of 
this equation is 

= . J, ... No (r - t) (17) 



■Q{E,t) 



(16) 



\dE/dt\' 



where 



dE 



\dE/dt\ 



(18) 



Figure 4. Spectrum of secondary positrons as a result of several 
successive capture events. The strength of magnetic field H = 3 
mG. 



and dE/dt is the rate of synchrotron or Coulomb losses. 

In the high energy range, where the synchrotron losses 
are essential, the spectrum is 



iV 



(1 - pEty 



E 



-l3Et 



e{l - l3Et) 



(19) 



E., 



(20) 



We see that independent of the injection spectrum the max- 
imum energy of positrons long after the injection is 

1 

- Jt 

where t is time after the injection. 

If the positrons are secondary, then their injection spec- 
trum has a cut-off in the low energy range at £cut — 30 MeV 
because of the threshold oi p — p reaction, and the energy 
distribution of positrons leaving the region of strong mag- 
netic fields looks like a succession of separated bunches (see 
FiglJJ. Here and below we take the injection spectrum of 
secondary positrons in the form 



Nq{E) = AtE'^e(E - 30 MeV) 



(21) 



where At ~ 5.4-10'^* MeV^ gives the total number of injected 
positrons ~ 3 • 10^^, wh ere the spectral index e quals —3 as 
in the Galactic disk (see iBerezinskii et al.|[l990l ). 
The bunch width is 



The time duration of the bunch (characteristic time during 
which the bunch crosses any energy E under the influence of 
synchrotron losses) is independent of E and for the spectrum 
(PT|) equals 

1 



Ti 



f3Ec 



(23) 



that gives = 3.6-10*yr for Ecut = 30 MeV and = 3 mG. 
Subsequent bunches are separated from each other by the 
periods of Ti ~ 10^ yr, the average periods of star capture. 

After the capture time equaled 10^ yr the maximum 
energy of positrons in the bunch shifts from oo to ~ 10 
MeV (see Eq. (|20p . and there are no positrons in the en- 
ergy range i5 > 10 MeV while the energy range below 10 



MeV contains permanently a number of bunches moving to 
the thermal region (see Fig. [4]). Then we conclude that the 
spectrum of positrons in the range above 10 MeV is strongly 
non- stationary (shown by the dashed line in Fig. |4]), as well 
as the emission generated by these particles. In contrary, the 
spectrum of positrons with energies < 10 MeV (shown by 
the solid lines in Fig. |4]) and the radiation produced by these 
positrons is quasi-stationary. 

When the rate of synchrotron losses drops down with 
positron energy, the following evolution of positrons in the 
range E < 1 MeV occurs under the influ ence of Coulomb 
losses. The rate of Coulomb losses is (|Havakawal Il964l : 
lGinzburel[l989l ) 



dE_ 
H 



2ne n 

'mcf3{E) 



In A 



(24) 



where log A is Coulomb logarithm, and /3(iJ) = v/c. For 
lorenz-factor 7 — — ^ -f 1 in a neutral medium 



log A ~ log [(7 - 1)(72 - 1)] + 20.5 , (25) 
while in completely ionized plasma 

logA~ log [7/ri] -1- 73.6 . (26) 



6 EMISSION PRODUCED BY 
NON-THERMAL POSITRONS 

Relativistic positrons with energies higher than 10 MeV gen- 
erate radio emission due to synchrotron losses and fluxes of 
in-flight annihilation and bremsstrahlung emission in the en- 
ergy range above 10 MeV. 

From Eqs. (fTi)l . fl^ and pTj) we calculated the time 
variations of radio emission from the GC at the frequency 
330 MHz. The result of these calculations is shown in Fig. 
[5^. As it was expected at the very beginning, just after the 
injection of high energy positrons, the flux of radio emission 
is as high as ~ 3 ■ 10^ Jy, but its value decreases very 

rapidly with time. The present time corresponds to periods 
when the radio flux drops down below the level about 1 kjy 
(dotted line), as the observations require. 
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Figure 5. a) Time- variations of the radio flux from the GC. Here 
^rnax° = 3 ■ lO'' Jy. b) Annihilation flux variations from the GC 
in the framework of non-stationary model. Here I^^^ = 8 ■ 10"'* 
photons cm~^s^^. 



Table 2. Flux of in-flight annihilation in case of non-stationary 
annihilation (photons cm~^s~^). 



Energy 


COMPTEL 2(7 


Ionized 


Neutral 


range 








(MeV) 








1-3 


2 X 10"* 


1.2 X 10""' 


2.6 X 10"^ 


3-10 


9 X 10-^ 


3.2 X 10"" 


3.6 X 10"'' 


10-30 


3 X 10~5 


3.7 X lO"'^ 


3.6 X lO"'^ 



Table 3. Flux of in-flight annihilation in case of stationary anni- 
hilation (photons cm~^s~-'^). 



Energy 

range 

(MeV) 


COMPTEL 2(7 


Ionized 


Neutral 


1-3 


2 X lO--* 


4.6 X 10-5 


1 X 10-* 


3-10 


9 X 10-5 


1.2 X 10-5 


1.4 X 10-5 


10-30 


3 X 10-5 


1.5 X 10-*^ 


1.4 X 10-^ 




1 2 3 4 5 6 7 

t(s) 



9 10 
X 10^^ 



Figure 6. Variations of the in-flight 10-30 MeV flux. The value 
of Fmax is presented in Table [2] 



The in-flight annihilation in the range 10-30 MeV is pro- 
duced by positrons whose intensity is non-stationary. There- 
fore, this flux from the GC is time-variable (see Fig. [6]) 

On the contrary, the flux of in-flight annihilation in the 
ranges 1-3 and 3-10 MeV is almost stationary as expected. 

The in-flight fluxes calculated for the cases of neutral 
and ionized medium and the COMPTEL 2a levels are pre- 
sented in Table (2] The flux of the in-flight annihilation in 
the range 10-30 MeV presented in the table corresponds to 
the moment when the radio at 330 MHz drops to 1 kjy. 

The bremsstrahlung flux in the energy range f rom 1 
to 30 MeV was calculated for the cross-section from iHaud 
l| 19971 ). Its values ranges from 5 x 10 ^ in the range 1 — 3 
Mev to 3 X lO-'^ at 3- 10 MeV and 1 x lO-'^ at 10-30 MeV. 
The flux units are the same as in Table (2] One can see that 
bremsstrahlung doesn't contribute much since lorentz-factor 
is rather small. 

From these tables one can see that in all energy 
ranges the calculated in-flight flux is below the COMP- 
T EL 2(7 level. Theref o re, in cont rast to the conclu sion 
of iBeacom and Yiiksell l|2006l ) and ISizun et al.l (|2006l ). we 



find that there is no problem with the injection energy of 
positron in the model with strong magnetic field at GC. Cor- 
responding combined spectra of in-fiight annihilation and 
bremsstrahlung emission for magnetic fields in the range 
0.1 — 3 mC are shown in Fig. [T] 



7 ANNIHILATION EMISSION OF THERMAL 
POSITRONS 

The characteristic time of annihilation for thermal positrons 
can be defined from 



Tann — 



E 



2y '^i^V^^^H-i)'^'^(^)^ 

Et 



(27) 
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E (IVleV) 



Figure 7. Combined in-flight annihilation and bremsstrahlung 
spectra for different values of magnetic fleld. COMPTEL limits 
are shown by heavy solid lines. 
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which is a function of the probability for thermal positrons to 
be in the energy range where the annihilation cross-section 
is high. Here temperature O ~ 1 eV as it usually assumed. 
For the gas density n > 1 cm""^ the annihilatio n time is 
about < 2.6 • 10^ yr (see iGuessoum et al.l |2005| ) which is 
smaller than the average time of star capture Ti and the 
bunch time Ti,. Then the process of annihilation emission 
is non-stationary since all positrons from a bunch annihi- 
late during the time shorter than the period between two 
neighbour capture events. 

The nondimensional equation for positrons has the form 
(here and below f{p) is the density of positrons with mo- 
mentum p) 



dt dp 



f 



op 



+ 



+novaiff + novacef = Q{p, t) 



(28) 



where dp/dt\^ and Dc are the rate of Coulomb losses and 
the momentum diffusion due to Coulomb collisions. These 
two terms form the Maxwellian spectrum of thermal par- 
ticles. The cross-sections aif and ace denote the in-flight 
annihilation of fast positrons and the charge-exchange anni- 
hilation process of thermal particles. 

The fluctuations of the annihilation flux expected in 
this case are shown in Fig. (Sja. The maximum intensity of 
the annihilation emission corresponds to the value of I^ax ~ 
8- 10~^ photons cm~^s~^. From this flgure one can see that 
there are long periods of time when the annihilation flux 
reaches its maximum value while the radio emission is below 
the 1 kjy level as the observations require, ft means that 
according to this model we see the annihilation flux from 
the bulge at its peak value. 

If the gas density no < 0.25 cm~^ then the annihila- 
tion time Tann is larger than Ti and the situation is quasi- 
stationary as follows from Eq. (|27p . Since the positrons from 
a single bunch do not annihilate completely for the period 
between two neighbour injections, then these positrons are 
accumulated in the thermal pool up to the level which pro- 
vides a quasi-stationary rate of annihilation. 

From Eq. ((28} we calculated the flux of annihilation in 
the quasi-stationary case which is shown in Fig. |S] by the 
solid line. We see that fluctuations of the annihilation flux 
is much smaller than in the non-stationary case (see Fig. 

m- 

The expected quasi stationary fluxes are presented in 
Table [S] The bremsstrahlung flux in this case is: 2 x 10~® 
at 1-3 MeV, 1 x 10"® at 3-10 MeV, and 3 x lO"'' at 10-30 
MeV. Again fluxes are below COMPTEL 2cr. 

An important conclusion from our calculations is that 
the only restriction of the model is the rate of synchrotron 
losses which should cool down positrons up to the energy 
< 10 MeV for the time < 10^ yr, see Fig. [1] The regions of 
permitted values of positron injection energies derived from 
the COMPTEL restriction for different values of the mag- 
netic field strength are shown in Fig. (9] the thick dashed line 
- ionized medium and the thick solid line - neutral medium. 
These region is below these lines. As one can see there are no 
restrictions for the injection energy from gamma-ray data if 
the magnetic field strength is high enough. 

However, as we discussed already in the case of high 
magnetic field we have an additional restriction - the fiux of 




Figure 8. Variations of the annihilation flux from the GC in 
the framework of quasi-stationary model: solid Une - without the 
effect of streaming instabiUty, dashed line - with the effect of 
streaming instability. Istat is the average value of the annihilation 
fiux. 




Figure 9. The multi-band limitations on the injection energy 
of positrons. Thin lines are the same as in Fig. \3\ Thick dashed 
line represents limitations from COMPTEL energy range in case 
ionized medium. Thick solid lines represents limitations from 
COMPTEL energy range in case of neutral medium. 



radio emission at 330 MHz. The range of energy values which 
are derived from the radio data are shown by thin lines. In 
that part Fig. [9] is the same as Fig. [3] each line correspond 
to different value of period of injections (in ascending order 
T = 1 X 10*, 5 X 10*, 10^ and 2 x lO'^ years). Thin solid 
lines correspond to injection in form of delta-function while 
thin dashed line correspond to power-law injection. Again 
as we mentioned already for strong magnetic field if period 
between two captures is long enough the injection energy 
can be very high. 

So combining the gamma and radio data we con- 
clude that the injection energy of positrons can be high if 
the model satisfies the tow conditions: the magnetic field 
strength is larger than 0.4 mG and the average period be- 
tween two neighbour captures is longer than 10^ yr. 
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Figure 10. Transformation of the bunch spectrum by the 
streaming instability. 



8 LANDAU DAMPING: QUASI-STATIONARY 
MODEL 



As it follows from Eq. (|19|) when a bunch of MeV positrons 
escapes the strong magnetic field region into the surrounding 
medium the positron distribution looks like a two-peak func- 
tion with one maximum at thermal energies and the other at 
the energies of the bunch (see Fig. [TO}. From the quasi- linear 
the ory it follows that s uch a particle distribution is unsta- 
ble llArtsimovich and Sagdeev 1979i . Ch.I, §1.16, p. 104, see 
also iLifshitz and Pitaevskiil Il98ll . Ch.III. §30, p. 243), and 
the flux of fast particles excites effectively plasma waves due 
the streaming instability . As a result a bunch of fast par- 
ticles with the energy E = Efj^nch is transformed into "a 
plateau" distribution in the energy region E < Etunch- This 
process of the bunch smearing occurs due to resonant wave 
excitation {v = ui^/k, where v if the particle velocity and k 
is the wave number). This process is described as diffusion in 
the velocity space with the diffusion coefficient Dst equaled 



Dst ~ 47r — -ntm 



(29) 



where nt and Vmax is the density and the velocity of bunch 
particles, Vmin is about the velocity of thermal particles 
and V is the current velocity of particles. The Langmuir fre- 
quency uip is 



•\/ 47rnoe^/m 



(30) 



where no is the density of backg round plasma. As numeri- 
cal c alculations showed (see, e.g. lArtsimovich and SagdeevI 
1 19791 . Ch.I, §1.16, p. 107) the bunch smearing was reached for 
several Langmuir times that in comparison with the model 
characteristic times is almost instantly. Calculations show 
that about one third of the bunch energy is transformed 
into the energy of excited plasma waves. 

The nondimensional equation for nonthermal positrons 
has the form in this case 



df d I dp 



The expected variations of annihilation flux from the 
GC are shown in Fig. [8] by the dashed line. On can see that 
effect of streaming instability does not change the flux value 
significantly. 



9 SPATIAL NON-UNIFORM MODEL 

We presented above our analysis of integrated fluxes from 
the GC. As we noticed mG magnetic fields occupy the inner 
sphere with the radius = 100 pc (i.e. its angular radius 
is < 1°) while the FWHM of the annihilation emission is 
about 6°. It means that in most of their lifetime positrons 
spend outside the central magnetic sphere. On the other 
hand, it follows from our analysis that positrons should be 
cooled down by strong magnetic fields up to the energy < 10 
MeV. These circumstances give restrictions for processes of 
positron propagation in the GC. 

Spatial variation of the positron distribution function 
with propagation te rms is described by the equation (see 
iBerezinskn et al.lll99dl '): 



J.2 



ON d 

^"aT ^ 'Be 



dE 
'dt 



N 



+ navN = QiE,r) , (32) 



Here dE/dt is the rate of total energy losses, and D, r is a 
spatial diffusion coefficient. We use a simple approximation 
for the spatial part of the injection function Q{r), assum- 
ing that positron sources are uniformly distributed inside 
the sphere of the radius ro, Q{r,E) — Q{E)6{ro — r). The 
function Q{E) is given by Eq. (|2H) but with the coefficient 
At — 5 ■ lO**® MeV'^/s as a result of averaging of injection 
processes over time (quasi-stationary injection). 
The boundary condition are 



(V. ■ iV)_o = 0, 
and 



- 



(Vp ■ iV), 







(33) 



(34) 



The loss term can be presented as 



dE 
It 



= (dE/dt)^^,^, + {dE/dt),c + 



+ {dE/dt)^^^^^ X e{r - T„) + 
+ (d^/rfi).«„s X 



(35) 



+novaiff + novacef = Q{p, t) 



(31) 



where (dE/dt) ^^^gj^ is the rate of synchrotron losses in the 
outer sphere r > th where the magnetic field strength equals 
H ~ 10~^ G, and (dE / dt) ^^^^ is the rate of energy losses 
inside the central sphere r < th losses in strong magnetic 
field // = 3 mC. 

The observed 6° FWHM of spatial distribution of an- 
nihilation emission corresponds to a sphere with the radius 
about 400 pc that requires the spatial diffusion coefficient be 
about Drr ~ 10^**cm^s~^. We notice here that propagation 
of MeV positrons in the GC is still very questionable (see 
e.g. lJean et~alll2009l ). 

Fast cooling by synchrotron losses is possible only in- 
side the sphere of strong magnetic field r < rn- Therefore, 
the radius of source region should be quite small, ro < rn- 
Only in this case positrons are cooled to MeV energies by 
synchrotron losses and then fill the 400-pc sphere. In Fig. [TT] 
we presented the calculated values of in-flight fluxes in the 
three COMPTEL energy ranges as a function of the source 
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Figure 11. Variations of the in-flight flux as a function of the 
source region radius. Here Fc is COMPTEL's 2a level. 



radius ro. These fluxes are normalized to the COMPTEL 2a 
level presented in Tabled As one can see extended sources 
with rorn are unable to satisfy the COMPTEL restrictions 
if the injection energy is 30 MeV. 



10 CONCLUSION 

We showed that: 



• UnUke 1 eacom and Yiiksell (|2006D and 1 izun et al.l 
l|2006l l. who restricted the injection energy of annihilating 
positrons at the value 1 — 7 MeV, we show that under the 
condition of strong magnetic fields H > 0.4 mG in the 
Galactic center this value can be much larger. This extends 
significantly a class of models acceptable for explanation of 
the annihilation emission from the GC. However, the nec- 
essary condition in this case is that the positrons should 
be cooled down by the synchrotron losses up to the energy 
< 10 MeV. Otherwise the expected radio flux and in-flight 
MeV emission from the GC exceed the level required from 
observations; 

• The main restriction of this model follows from the 
value of observed radio flux from the GC which equals 1 
kjy at the frequency 330 MHz. In order to satisfy this con- 
dition the peak of annihilation emission should be shifted in 
time to the moment when the radio flux falls down below 
the observed level; 

• The energy spectrum of annihilating positrons looks 
like a number of bunches produced by subsequent capture 
events. In the case of the secondary origin of positrons the 
structure of the bunch spectrum is almost quasi-stationary 
in the energy range < 10 MeV, and strongly non-stationary 
in the range above 10 MeV; 

• On the other hand, the expected flux of the annihilation 
emission is quasi-stationary if the density of the background 
gas is < 0.25 cm~^. 

• Time variations of the in-flight annihilation in the en- 
ergy ranges 1-3, 3-10, and 10-30 MeV are quite small, and 
the important point is, their values are smaller than the 
COMPTEL 2a level; 

• Our calculations show that the size of the source region 
should be smaller than the radius of the sphere filled with 



strong magnetic field. Otherwise the in-fiight fiux exceeds 
the COMPTEL 2a level. 

It would be interesting to find traces of similar processes 
from other galaxies. If positrons injecti on is related to stel - 
lar capture events as it w as suggested bv lCheng et al.l (j2006l ) 
and lCheng et al.1 |20o3) one can expect to find some traces 
of these processes in galactic nuclei with recent capture 
events in the form of high radio flux. However, galactic 
nuclei with high X-ray flux that indicates active accretion 
processes there do not show such a powerful radio emission 
(jWong et al.ll20o"7l ). This is not surprising since the peak of 
radio emission is expected in our model some time after the 
capture. The delay time is about 10^ years as follows from 
the characteristic time of p — p collisions which produce high 
energy electrons: 



{nap 



10' yr 



(36) 



for the gas density of molecular clouds n ~ 10'* cm ^ . 

The flux of X-ray emis sion from the accretion disk scales 
with time like (|Reeslil988l ') 



( — ) 



-5/3 



(37) 



where tmin — 0.2 yr. So by the time radio emission peak 
reaches the X-ray emission decreases by a factor of 10^ from 
its initial value ~ lO'** erg s~^. 

Thus, we expected a radio flux with the value ~ 10* Jy 
when the X-ray flux is already unseen. Galactic nuclei with 
such a brightness are c ommonly presented in the Universe 
iHveop Lee et al.ll2009l ) and what is interesting that in some 
cases there is no activity in other energy ranges. That falls 
well into our model though these speculations cannot be 
considered as direct evidences in favor of our model. 
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